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We report quantum-mechanical close-coupling calculations for rotational state resolved 
cross sections for P-H2+HD collisions. The low temperature limit of P-H2+HD is investi- 
gated, which is of significant astrophysical interest in regard to the cooling of primordial 
gas and the interstellar media. Sharp resonances have been reproduced in the cross sec- 
tions of some transition states at very low kinetic energies, Ekin ~ 10"^ eV. 

1. Introduction 

Energy transfer collisions between H2 and HD molecules, where H is hydrogen and D 
is deuterium, is of fundamental importance for the astrochemistry of the early Universe 
and the interstellar medium [[Il[2l[3lll[5l[6l[71[8l[9l[T0l[ni[T2l[T3]. HD represents the 
second most abundant primordial molecule after H2 and plays a significant role in the 
cooling of the primordial gas. The abundance of deuterium is low relative to hydrogen 
(D/H?a 10-^), however the HD/H2 ratio can be enhanced by an approximate factor of 
10^ due to chemical fractionation [ IH E] • 

It has been shown, that in the framework of the standard cosmological model the radi- 
ation temperature is higher than the matter temperature and molecules become a heating 
source for the gas. At higher temperatures H2 molecules dominate the heating, however 
HD molecules dominate the kinetic process at lower temperatures. In addition, the HD 
molecule is especially important due to its permanent dipole moment and lower rota- 
tional constant, which makes the molecule to be an efficient coolant at low temperatures: 
T < lOOK [[l3]. When H2 molecules are inefficient HD become important in cooling the 
primordial gas and the interstellar medium. This ability of HD makes these molecules to 
be very attractive. 

Therefore knowledge of the ro-vibrational excitation and de-excitation thermal rate 
constants in the molecular H2+HD low energy collisions is of fundamental importance 
in understanding and modeling the energy balance within primordial gas. However, to 
accurately model the thermal balance and kinetics of such systems one needs accurate 
state-to-state cross-sections and rate constants k^j^iji{T). 

Experimental measurement of quantum state resolved cross sections and rate coeffi- 
cients is a very difficult technical problem. On the other hand, accurate theoretical data 
require precise potential energy surfaces and reliable dynamical treatment of the collision 



*rasultanov@stcloudstate.edu; r.sultanov2@yahoo.coni 
'''dcguster@stcloudstate.edu 



2 



R.A. Sultanov and D. Guster 



processes. The first attempt to construct a realistic full-dimensional ab initio PES for 
the H2-H2 system was done in works [ HH [T5] , and the potential was widely used in a 
variety of methods and computation techniques. Currently the deuterium chemistry of 
the eraly Universe has been extensively studied by many researchers. Here we would like 
to mention works [ [71 [8] and [ |3] , where the authors carried out quantum- mechanical 
state-resolved calculations for low and very low temperatures down to 10 K [|3], however 
using of the old and modifyed H2-H2 surface [E]. 

Nontheless the importance of the H2+HD in the astrophysical problems makes it vital 
to carry out new calculations for the system with recently published PES [US]. This is 
done here in. Moreover, with the current calculations we extend our recent experience for 
H2+H2 at low and very low kinetic energies [[T71 [18] . 

As was mentioned above a new extensive study of the hydrogen-hydrogen surface has 
been recently reported by Boothroyd et al. [[IS], where the potential energies have been 
represented at 48180 geometries respectively with a large basis set at the multireference 
configuration interaction level. In this work we provide the first calculations describing 
collisions of rotationally excited H2 and HD molecules using the BMKP PES 



The scattering cross sections are calculated using a non reactive quantum-mechanical 
close-coupling approach. In the next section we will briefly outline the method and present 
results for the cross sections for rotational de-excitation of HD in low energy collisions 
with p-B.2. We compare our results with some previous investigations. Conclusions are 
presented in Section 3. 

2. Method and Results 

In this section we provide a brief outline of the quantum-mechanical close-coupling 
approach used in our calculations. The basis for this methodology was developed in work 
[[19]. The HD and H2 molecules are treated as linear rigid rotors. For the considered range 
of kinetic energies of astrophysical interest and for practical astrophysical estimations the 
rotor model is considered to be adequate for P-H2-I-HD collisions [IB]. 

The 4-atomic H2-HD system is shown in Fig. 1. It can be described by six independent 
variables: Ri and R2 are interatomic distances in each hydrogen molecule, Bi and B2 
are polar angles, $2 is torsional angle and R3 is intermolecule distance. Let us introduce 
M12 = (mi+m2)(m3 + m4)/(mi + m2 + m3 + m4) and /ii(2) = mi(3)"^2{4)/("^i{3) + "^2(4)), 
where the first one is a reduced mass of the pair of two- atomic molecules 12 and 34 and the 
second ones are reduced masses in hydrogen molecules. As we mentioned, the hydrogen 
molecules are treated as linear rigid rotors, that is distances Ri = 0.7631 a.u. in H2 and 
R2 = 0.7668 a.u. in HD are fixed in this model. We provide a numerical solution for the 
Schrodinger equation for a 12 + 34 collision in the center of the mass frame. 

The cross sections for rotational excitation and relaxation phenomena can be obtained 
directly from the S-matrix. In particular, the cross sections for excitation from jij2 — > ^2 
summed over the final m'im'2 and averaged over the initial mim2 corresponding projections 
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of the H2 molecules angular momenta ji and j2 are given by 



-(A^&nn,.)- E (2j + i)|&„,-si,(£)P. (2) 

The kinetic energy is e = E — Biji (ji + 1) — -62^2 ( j2 + 1) • Here E is the total energy in the 
system, Bi = 60.8 cm~^ and B2 = 44.7 cm~^ are the rotation constants of the colliding 
HD and H2 molecules respectively, J is total angular momenta of the 4-atomic system, 
a = (jij2ji2^), where ji + j2 = ji2 and ji2 + L = J, k^a' = ^M^iE + Ea- E^if/^ is the 
channel wavenumber and Ea{a') are rotational channel energies. 

We apply the hybrid modified log-derivative-Airy propagator in the general purpose 
scattering program MOLSCAT [ [20] to solve a set of coupled second order differential 
equations for the unknown radial functions U^^{K) 



cP L{L + 1) 



dR^ R^ 

\V{n,r2,RM'J'{ri,r2,R) > UiF {R)dhdf2dR, (3) 

We have tested other propagator schemes included in the MOLSCAT code. It was found, 
that other propagators can also produce quite stable results. 

The log-derivative matrix is propagated to large i?-intermolecular distances, since all 
experimentally observable quantum information about the collision is contained in the 
asymptotic behaviour of functions U^^\R 00). The numerical results are matched to 
the known asymptotic solution to derive the physical scattering ^-matrix 



where kaa' = 2Mi2{E + Ea — Ea'Y^'^ is the channel wavenumber, Ea(a') are rotational 
channel energies and E is the total energy in the 1234 system. The method was used for 
each partial wave until a converged cross section was obtained. It was verified that the 
results are converged with respect to the number of partial waves as well as the matching 
radius, Rmax, for all channels included in our calculations. 

The new BMKP PES [IS], which is used in these calculations, is a global six- dimensional 
potential energy surface for two hydrogen molecules. It was especially constructed to rep- 
resent the whole interaction region of the chemical reaction dynamics of the four-atomic 
system and to provide an accurate as possible van der Waals well. In the six-dimensional 
conformation space of the four atomic system the conical intersection forms a complicated 
three-dimensional hypersurface. 

Because the BMKP PES [ [16] uses cartesian coordinates to compute the distances 
between four atoms, we needed to devise a fortran program, which converts spherical 
coordinates used in the close coupling method [ 20j to the corresponding cartesian coor- 
dinates and computes the distances between the four atoms followed by calculations of 
interatomic interaction forces. 

The four atomic system is shown in Fig. 1. Let us introduce the Jacobi coordinates 
{Ri, i?2, -R3} and the radius- vectors of all four atoms in the space- fixed coordinate system 
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OXYZ: {ri, r2, ^3, ^4}- We apply the following procedure: i?3 is directed along OZ axis, 
the center of mass of the HD molecule is brought into coincidence with the center of 
OXYZ, and the top of the i?3 is directed to center of mass of H2, as shown in Fig. 
1. Now it is apparent, that i?3 = {i?3, 83 = 0, <l'3 = 0}, -Ri = ri — r2, R2 = — rs, 
and fi = aRi and ^2 = (1 — a)Ri, where a = m2/{mi + m2). Next, without the loss of 
generality, we can adopt the OXYZ system in such a way, that the HD interatomic vector 
Ri lies on the XOZ plane. Then the angle variables of Ri and R2 are: Ri = {Qi, $1 = vr} 
and i?2 = {©2, ^2} respectively. 

Now one can see, that the cartesian coordinates of the atoms of the HD molecule are: 
ri = {xi = aRismQi,yi = 0,zi = a-Ricos0i}, r2 = {^2 = (1 — q;)-Ri sin 9i, ?/2 = 
0,^2 = —(1 ~ a)-RiCos0i}, and in turn for the H2 molecule we have: = {x^ = 
-(1 - (3)R2 sin 62 cos $, ?/3 = -(1 - /5)i?2 sin 02 sin $, = Rs - {1 - /3)i?2 cos 62} and 
r4 = {x4 = /3-R2 sin 62 cos $, 1/4 = /3i?2 sin 62 sin $, 2:4 = R3 + PR2COSQ2}, because = 
^3 — (1 — P)R2 and r4 = ^3 + PR2, where /3 = m^/^m^ + 1714). In such a manner the 
cartesian and the Jacobi coordinates are represented together for the four-atomic system 
H2+HD. 

Before our production calculations we carried out a large number of test calculations 
to insure the convergence of the results with respect to all parameters that enter into the 
propagation of the Schrodinger equation. The same calculations were also done in our 
previous works [ITTJIIH] for the 0-/P-H2+H2 collisions, which involved the intermolecular 
distance R, the total angular momentum J of the four atomic system, the number of 
rotational levels to be included in the close coupling expansion Niyi, and others (see the 
MOLSCAT manual 

We reached convergence for the integral cross sections, cr{j[, j2, jij2, e), in all considered 
collisions. For example, for R we used from Rmin = 1 A to Rmax = 22 A, we also apphed 
a few different propagators included in the MOLSCAT program. We obtained convergent 
results in all cases. 

Below we present our calculations for the state-resolved cross sections crOijj.) iiJ2)(e) 
in the collision (1). Fig. 2 shows our data for the (ji = 1, j2 = 0) — {j[ = i'2 = 0) 
and (ji = 2,j2 = 0) ^ {j[ = l,j2 = 0) quantum transitions together with the results of 
Schaefer's work [[3], which applied a modified version of Schwenke's potential [HI] and 
different dynamical quantum-mechanical approach. As can be seen for these transition 
states the agreement between the two different calculations is excellent. We reproduced 
sharp resonances in the low velocity region, which are very important in cooling of the 
astrophysical media. In addition to Fig. 2 we separately show the cross section for only 
(ji = l,j2 = 0) — s> {j[ = i2 = 0) transition states in Fig. 3. This graph can be directly 
compared with the corresponding cross section from the Schaefer work [[3], please refer 
to figure 6 of that paper [ [3] . One can see, that even small peculiarities at around 250 
m/s and 450 m/s of the cross section behaviour are reproduced in details. 

In the Fig. 4 we present the energy dependence of the cross sections for the (ji = 2, 72 = 
0) {j[ = i'2 = 0) and (ji = 1,^2 = 2) ^ (j[ = 0, j2 = 2) quantum-mechanical transition 
states. For the last transition we obtain very good agreement with the corresponding 
results of work [ |3]. However, in the upper graph of Fig. 4 one can see rather large 
differences in the low velocity region. Our calculations revealed, that the disagreement 
ranges up to about 100 %. 



Rotational energy transfer in p-B.2 +HD 



5 



Finally, significant differences are found in the cross sections from the higher values of 
the transition states. Fig. 5 presents our results again together with Schaefer's data from 
m for (ji = l,j2 = 2) ^ ij[ = 2,j^ = 0), (ji = l,j2 = 2) ^ ij[ = = 0) and 
(ji = 1, j2 = 2) ^ {j[ = j'2 = 0). In this case the disagreement ranges up to one order of 
magnitude. 

Also, as can be seen from the graphs we calculated few new resonances for each of these 
transitions at a velocity of about 1200 m/s. The biggest value of these resonances is in the 
cross section (ji = l,j2 = 2) — >■ {j[ = 2,j2 = 0). We show this separately in more detail 
in Fig. 6. The value of the resonance is relatively large and it also might be applicable in 
important astrophysical processes, such as coUisional cooling. 

Through this analysis we can now conclude, that the new global BMKP PES can 
reproduce general behaviour of all considered cross sections in the P-H2+HD collision. 
For the lower quantum states we obtained sufficient agreement with previous calculations. 
However, for transition states from higher values, for example, (ji = 2, j2 = 0) the BMKP 
PES provides rather low cross sections relatively to the corresponding results of work [|3]. 

3. Conclusion 

In this letter the state-to-state close-coupling quantum-mechanical calculations for rota- 
tional excitation and deexcitation cross sections of the P-II2+IID collision are presented. 
The linear rigid rotor model for the H2 and HD molecules is applied. The global and 
newest BMKP surface for the H2-H2 system has been appropriately adopted for the cur- 
rent P-II2+IID collisions. A test of convergence and the results for cross sections with the 
BMKP PES are obtained for a wide range of kinetic velocities including very low values 
down to 10 m/s. 

Our results revealed, that for low quantum transition states the BMKP surface provides 
cross sections very close to those obtained in previous works [ El C] , where the authors 
adopted Schwenke's old H2-H2 PES [Ej. However, for some higher quantum states we 
found significant disagreements with previous results [|3] (see Fig. 5). Additionally, in our 
calculations some new resonances are found in the 1300±100 m/s region for transition 
states from j2 = 2 and ji = 1. It was found, that for the specific transition (^2 = 2, ji = 
1) ~^ = 2, ii = 0) the value of the resonance is relatively large in the cross section, and 
it may even stronger infiuence, for example, the cooling processes in primodial gas and 
interstellar media. Further detailed calculations for higher quantum transition state cross 
sections and corresponding thermal rate constants kj^j^^ji^ji^{T) for 0-/P-H2+HD collisions 
at low and very low kinetic energies are in progress in our group. 

In conclusion, we would like to point out here, that the ultralow energy sector (T < 
1 /iK) is of crucial importance now in connection with the recently achieved molecular 
Bose-Einstein condensates [ |2T1 [22l [23l [2ll [251 126]. Therefore, it should also be inter- 
esting and useful to apply the current time-independent, quantum- mechanical approach 
to investigate isotope effects in molecular hydrogen collisions at such ultralow energies. 
Further, the methodology could be used to carry out new calculations for the important 
ultralow collisions considered in the works [ [2Z1 ESI EH] , as well as check the conclusions of 
the recent work [ [30] involving scattering properties of weakly bound dimers of fermionic 
atoms. 
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Figure 1. Four-atomic coordinates for the p-H2(j2)+HD(ji) collision used in this work. 
Ri = 0.7631 a.u. and R2 — 0.7668 a.u. are fixed interatomic distances in each hydrogen 
molecule HD and H2 respectively, ©i and ©2 are polar angles of vectors Ri and R2 
respectively, $2 is torsional angle and R3 is the inter molecular vector, which connects 
the center of masses of the molecules. Vectors ri,r2,r3 and represent the cartesian 
coordinates of the four atoms in the space-fixed OXYZ coordinate system. 
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Figure 2. Rotational state resolved integral cross sections for p-H2(j2) + HD(ji) — >• 
^2(^2) +HD(j^). Initial states of HD and H2 molecules are ji = 1 and j2 = respectively 
and corresponding final states are j[ = J2 = 0. In the bottom plot: ji = 2, ^2 = and 
j[ = 1, ^2 = 0. Calculations are done with the BMKP PES (bold lines), triangles left are 
corresponding results from work [[3]. 




Figure 3. Details of the (ji = 1,^2 = 0) - 
cross section in p-B.2{j2 = 0) + HD(ji = 1) 
those corresponding cross sections from wo 



{j[ = 0,^2 = 0) rotational transition state 
. H2(j2 = 0) + HD(j; = 0) to compare with 
[[3] figure 6. 
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Figure 4. Rotational state resolved integral cross sections for p-H2(j2)+HD(ji) — > H2(j2) + 
HD(j^). Initial states of HD and H2 molecules are ji = 2, j2 = respectively and 
corresponding final states are j[ = ^2 = 0. In the bottom plot: ji = 1, j2 = 2 and 
j[ = 0, = 2. Calculations are done with the BMKP PES (bold lines), triangles up are 
corresponding results from work [ |3] . 
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Figure 5. Rotational state resolved integral cross sections for p-H2(j2)+HD(ji) H2(j2) + 
HD(j^). The initial states of HD and H2 molecules are ji = 1 and j2 = 2 respectively and 
corresponding final states are j[ = 2,^2 = in the upper, j[ = 1, j!^ = in the midle, and 
j'l = ^2 = the bottom plots. Calculations are done with the BMKP PES (bold lines), 
diamonds are corresponding results from work [ [3] . 
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Figure 6. Sharp resonance in the (ji = 1, j2 = 2) ^ {j[ — 2, = 0) rotational transition 
state cros section of p-H2(j2) + HD(ji). The kinetic velocity ranges from 1000 m/s to 
1800 m/s in this graph. 



